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Abstract. We study fundamental spectral properties of random 
block operators that are common in the physical modelling of meso- 
scopic disordered systems such as dirty superconductors. Our re- 
sults include ergodic properties, the location of the spectrum, ex- 
istence and regularity of the integrated density of states, as well as 
Lifshits tails. Special attention is paid to the peculiarities arising 
from the block structure such as the occurrence of a robust gap 
in the middle of the spectrum. Without randomness in the off- 
diagonal blocks the density of states typically exhibits an inverse 
square-root singularity at the edges of the gap. In the presence 
of randomness we establish a Wegner estimate that is valid at all 
energies. It implies that the singularities are smeared out by ran- 
domness, and the density of states is bounded. We also show Lif- 
shits tails at these band edges. Technically, one has to cope with 
a non-monotone dependence on the random couplings. 



I. Introduction 

Random block operators play an important role in the mathemat- 
ical modelling of superfluid fermions in a random environment and 
are thus relevant for mesoscopic disordered quantum systems such as 
dirty superconductor devices. They arise in the Bogoliubov-de Gennes 
equation 



B* -H ) ) \if>- 

that is, the eigenvalue problem for the quasi-particle (or excitation) 
states (J+) in a mean-field approximation of BCS theory |dG] . Without 
loss of generality we have assumed that the chemical potential equals 



zero in (1.1). The 'particle' and 'hole' components ip + and of the 
quasi-particle state belong to the single-particle Hilbert space %. The 
self-adjoint single-particle Hamiltonian H = H* and the so-called the 
pair potential or gap function B are linear operators on %. The overbar 



in (1.1) denotes complex conjugation. 



Following Altland and Zirnbauer |A1Z] one can classify all block op- 



erators that arise in (1.1) according to their behaviour with respect to 
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time-reversal and spin-rotation symmetry. In this paper we will focus 
on random block operators of the form 

H:=(f_*J (1-2) 



B —H 

with both H and B self-adjoint. This choice corresponds to symme- 
try class CI of [A1ZJ and describes physical systems for which both 
time-reversal and spin-rotation symmetry hold. Since the Bogoliubov- 
de Gennes equation results from a mean-field approximation, the ex- 
pressions for the operators H and B should be determined from self- 
consistency requirements. For disordered systems the discrete Ander- 
son model in d dimensions, 

H:=A + V on H = f(Z d ) } (1.3) 

is a generally accepted effective description for this, see e.g. |ViSFj . 
(Choosing H as a random matrix from a suitable ensemble would be 
another |AlZj .) Here, 

(A^)(j):= m (1.4) 

i&L d :\i-j\=\ 

for all j G Z d and all ip G £ 2 (7j d ) is the centred discrete Laplacian 
and the random potential V amounts to multiplication by independent 
and identically distributed, real- valued random variables {V(j)}j &Z d 
according to (Vip)(j) := V(j)if)(j). 

The form of the gap operator B, which should also be determined by 
self-consistency, depends on the pairing mechanism. For s-wave (a.k.a. 
conventional) superconductors B is a multiplication operator in posi- 
tion space £ 2 (7j d ). Homogeneous s-wave superconductors are described 
by a multiple of the identity operator, B = (3 1 with a self-consistently 
determined parameter j3 > 0. Disordered s-wave superconductors are 
often described by an effective random multiplication operator B = b 
[dOl IZl IViSF] . Here (bip)(j) := b(j)if>{j), where {b(j)} jeZ d are inde- 
pendent and identically distributed real-valued random variables. In 
addition, the 6(j)'s are often required to be independent of the V(j)'s. 
Our main results in Sections [5] and [6] will be proved in precisely this 
setting. 

Non-diagonal gap operators B occur in the modelling of cl^^-wave 
superconductors. For example, the momentum-dependent interaction 
of Cooper pairs leads to B = (3(A^ ] - A^) for homo geneous supercon- 
ductors in two dimensions |DL1 lAlSZj , where /3 > and denotes 
the one-dimensional centred discrete Laplacian in x-, resp. y-direction. 
Our results may also be of relevance for inhomogeneous (i-wave super- 
conductors [Zj IViSFj , if we use models with diagonal disorder for these 
materials. 

The plan of this paper is as follows. In Section [2] we study basic 
spectral features of block operators H that are of the general form 
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Figure 1 . Comparison of the density of states of a ran- 
dom block operator H without (left) and with (right) 
randomness in the off-diagonal blocks B. The singular- 
ities at the inner band edges are smoothed out by the 
randomness in B. 



(1.2). Among others we establish the existence of a robust spectral 



gap of H in Proposition |2.10[ We interpret the robustness of the gap 
in the context of Anderson's theorem |An| IBVZ] in Remark 2.11 (iii)| 

Section [3] briefly discusses an important special case of (1.2), namely 
constant off-diagonal blocks B = (31. This serves to expose a typical 
phenomenon in the absence of disorder: the density of states of H 
suffers from an inverse square-root singularity at the inner band edges 
in every dimension d G N. The singularity is robust in the sense that 
it always shows up unless the density of states of H vanishes at energy 
zero - the location of the chemical potential. 

We introduce the main objects of this paper, ergodic random block 
operators, in Section [4j The basic spectral consequences of ergodicity 
are also explored there. This includes the location of the almost-sure 
spectrum, as well as the definition and self-averaging of the integrated 
density of states of H. 

In Section [5] we show that the density of states exists for suitable 
ergodic random block operators and that it is bounded. This is the 



content of Theorem |5.1[ the main result of this paper, and follows 
from a Wegner estimate. Figure [T] compares this situation to the one 
with the singularity at the inner band edges for a constant B — (3 1 as 
in Section |3j Randomness in the off-diagonal blocks smooths out the 
singularities. We stress that the Wegner estimate of Theorem 5^ holds 
for a block random Schrodinger operator with a sign-indefinite single- 
site potential of mean zero. In contrast, for ordinary (i.e. non-block) 
random Schrodinger operators such a result is still missing despite a 
lot of recent efforts IVel 



Finally, Section [6] establishes Lifshits tails for the integrated density 
of states of H at the inner band edges. Theorem 6A embodies the 
second main result presented here. 
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2. Basic properties of block operators 

In this section we explore some fundamental properties of rather 
general self-adjoint block operators 

H:=(2 B „) (2.1) 

on the Hilbert space H 2 := W © %. We equip % 2 with the scalar 
product 

:=(il>i,<Pi) + (il>2,<P2), (2.2) 
where (•, •) stands for the Hilbert space scalar product of TL and \& : = 

(S> $ : = © G U2 - We write 1*11 := dl^ll 2 + ll^ll 2 ) 172 for the 
induced norm. 

Later we will be mainly interested in the case where the self-adjoint 
operator if is a discrete Schrodinger operator on TL = C 2 {'L d ) and B is 
a multiplication operator with some real-valued function b on 7L d . In 
this situation both operators H and B are frequently bounded, so that 
H is unambiguously well defined as a self-adjoint operator. However, 
we can treat unbounded operators as well, as one can read off from the 
following assertion. 

Proposition 2.1. (i) Let H,B be self- adjoint, assume that 
dom(B) fl dom(if) is a core for H and that domdH] 1 ^ 2 ) C dom(B). 
Then H is essentially self-adjoint on ( dom(_B) ndom(if)) © ( dom(B) D 
dom(H)). 

(ii) Let H be self-adjoint, let B be symmetric and H -bounded with 
bound strictly smaller than one. Then H is self-adjoint on dom(if ) © 
dom(if). 

(iii) Let B be self-adjoint, let H be symmetric and B-bounded with 
bound strictly smaller than one. Then H is self-adjoint on dom(£?) © 
dom(.B). 

Proof. The assertions are special cases of Thm. 2.6.6 and Prop. 2.3.6 
in p. □ 

Without further mentioning we will assume in the rest of this paper 



that at least one of the three situations described by Proposition 2A_ 
applies, thereby ensuring self-adjointness of H. Next we compile some 
basic structural properties of the spectrum of H. 

Lemma 2.2. The operators H and H' := ^ B J are unitary equiva- 
lent. 

Proof. Define the unitary involution Ui := r^(i_x)> then H' = 
UiHUi*. X ~ □ 

Lemma 2.3. The spectrum o/H is symmetric around ; i.e. 

spec(H) = -spec(H). (2.3) 
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In particular, if H\I/ = for some £eK and \l> = Oi,^ 2 ) T e U 2 , 
then = -E^ , where § = (^ 2 , -^i) T - 

Proof. We define the unitary transformation U2 := (_ 10 ) 011 
which obeys U^ 2 = -U 2 , and observe U 2 HU^ = — H. □ 

Remarks 2.4. (i) It follows that the spectrum of H 2 has multiplic- 
ity at least 2, except possibly at zero, and that 

spec(H) = {E ER\E 2 e spec(H 2 )}. (2.4) 

(ii) The anti-symmetry of H under the transformation U2 is known 
as 'particle-hole' symmetry. 

Lemma 2.5. The operator H 2 is given by 

2 ( H 2 +B 2 [H,B] 



and unitarily equivalent to © K + on H 2 , where 

K± := H 2 + B 2 ±i[H,B]. (2.6) 

Proof. This follows from an explicit computation and the observation 
that 

U 3 H 2 U;= ( K - ° ) (2.7) 



K , 

for the unitary U 3 := \ ( \ \ ) . □ 

Remark 2.6. Introducing the annihilation operator a := H — i B on %, 
we can write K + = aa* and K_ = a* a. Thus, the spectra of K + and 
K_ differ at most by {0}. In fact, we have 

( K ol) = C°of U 3 HU; = i(° -„»*). (2.8) 

A direct calculation also shows the next assertion. 

Lemma 2.7. Suppose there exists a unitary involution U = U* = U^ 1 
on H such that HU + UH = and [B, U] = 0. Then U := ^= ( J _^ 

zs unitary on H 2 and UHU* 



H+UB 
ii-C/B 



Later we will use the following particularisation of Lemma 2.7 



Corollary 2.8. Assume B = b is the maximal self-adjoint multiplica- 
tion operator by the function b : Z d — > R and A t/ie centred discrete 

Laplacian as defined in (1.4). T/ien £/ie operator H = ^ _ fe A J zs mhj- 
tarily equivalent to H + © where H± := A ± (— 1) J 6. 

Proof. Choose U := (— that is 

(WW) ■= (-l)V(j) := (-1)^HU) (2-9) 
for all j = (ji, . . . , jd) G Z d and ip G £ 2 (Z d ), and apply Lemma 2.7 □ 
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Remarks 2.9. (i) It is essential for the validity of Corollary 2.8 
that A contains no diagonal terms. 

(ii) Lemma 2.5 and Corollary 2.8 imply that H 2 is both unitarily 
equivalent to K- © K + and also to if 2 © Hi . We have UK + U = K_ 
and UH + U = —H— Note also that none of the operators H^_ or i? 2 
alone is unitarily equivalent to K + or K_. 



iii) From Corollary 2.8 we can infer some information on the lo- 



cation of the spectrum of discrete Schrodinger operators with periodic 
potentials. For example, if b is the constant function equal G R, 
then the corresponding operators H± are unitarily equivalent. Hence 
spec(H + ) = spec(i/_) = spec(H). Thus it will follow from Proposi- 
tion 2.10 below, that the interval ] — (3,(3[ is a spectral gap for the 
operator A + (— l)- 5 /3. 

Finally, we are concerned with locating the spectrum of general H 
and with the occurrence of spectral gaps that arise from the special 
block structure. 

Proposition 2.10. (i) spec(H) C [ - \\H\\ - \\B\\, \\H\\ 
(ii) Suppose there exists A ^ such that H ^ XI. Then 
spec(H) n ] - A, A[= 0. 
Suppose there exists (3^0 such that B ^ (31. Then 

spec(H) n ]-p,p[= 0- (2.11) 
Suppose there exist A, ^ such that H ^ XI and B ^ (31. 

spec(H) n ] - (A 2 + /3 2 ) 1/2 , (A 2 + /3 2 ) 1/2 [ = 0. (2.12) 



+ \\B 



(2.10) 



m 



Then 



Remarks 2.11. (i) The endpoints of the interval in Proposition 
2.10 (i)| are sharp upper and lower bounds for the maximum and mini- 
mum of spec(H). This can be seen by choosing B as the multiplication 
operator by the function b = (— l)- 7 in Corollary 2.8 



(ii) The statements of Proposition 2.10(H), |(iii)| and (iv)| remain 



true if one replaces H and B by — H and 

(iii) According to Proposition 
gap of size at least 2(3 



B in the assumptions. 

H has always a spectral 
no matter what H is like. This statement 



2.10 iii 



can be interpreted in the context of Anderson's theorem \An\ IBVZj : 
Anderson argued that adding non-magnetic impurities to an s-wave 
superconductor should have only little (or at best no) effect on the 
gap operator B. In view of Proposition 2.10 (iii), this implies stability 



of the gap under non-magnetic impurity doping and thus leads to the 
experimentally observed insensitivity of superconductivity in this case. 
As an aside we mention that the situation is totally different when 
doping d x 2_ y 2-w&ve superconductors by non-magnetic impurities. Here 
it is known that disorder leads to the breaking of Cooper pairs. 
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Proof of Proposition 2.10 . The statement in follows from Lemma 
2^|and H^ill «C (\\H\\ + ||.B||) 2 . 

To prove (iv) we define H := H — A ^ and B := B — (3 ^ 0. Then 



we get for the operators K± in Lemma 2.5 

K±=H 2 + B 2 ± i[H, B] + X 2 + 2XH + /3 2 + 2/3B ^ A 2 + /3 2 , (2.13) 
where we used that the linear terms in A and /3 are manifestly non- 



negative, as is the sum of the first three terms on account of Remark [276 
Hence, H 2 ^ A 2 + (3 2 . 

and [( 



The proofs of 
same lines. 



m 



are simpler and proceed along the very 

□ 



Part (ii) of Proposition 2.10 allows for a generalisation towards dif- 



ferent diagonal blocks. 

Lemma 2.12. Let H\ be a self-adjoint operator on H, such that the 

( H B 

exists A ^ such that both H ^ A and H\ ^ A. Then 



block operator Hi 



is self-adjoint on H 2 . Suppose there 



spec(Hi) n ] - A, A[= 0. 



(2.14) 



Proof. We define H := H - A ^ and H x := Hi - A ^ 0. An explicit 
computation gives 

T 2 _ ( H 2 +B 2 HB—BHi \ _ / X 2 +2XH N ' " " x2 
BH—HiB H 2 +B 2 ) ~ \ o X 2 +2\H 1 



H 2 



H B 

B -H 1 



> A 2 1 



(2.15) 
□ 



3. The case of constant diagonal B 

As a warm-up and in order to expose some typical features we con- 
sider first the simple case where B = (3 1 is a constant multiple of the 
identity. In particular, the operators H and B commute. As a conse- 
quence, the spectral theory of H can be reduced to the diagonalisation 
of a 2 x 2-matrix and the spectral theory of H. 

Proposition 3.1. Let H be a self-adjoint operator on H, let f3 G M\{0} 
and consider 

H ==(£-*) 

on H 2 . Then 

(i) spec(H) = {E± := ±(E 2 + /3 2 ) 1 / 2 : E G spec(H)}. 

(ii) If Hip = Eip for some E G spec(if) and (p G H., then H$± = 
E±&± for the non-normalised vector Q± := \a-\tE -E) ) ^ ^ ' 
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Figure 2. Comparison of the density of states of a ran- 
dom block operator H with B = (left) and B = 01 
(right). 



(iii) Fix also E e R. Then there exists a constant C = Ce,p > 
such that if \\(H — E)cp\\ ^ e for some <p E % and e > 0, then 

\l(n-E±)$ ± \l^Ce. (3.2) 



Proof. Part (i) is a corollary of Lemma 2.5, Part |(ii) follows from an 
explicit computation. As to Part |(iii) we observe 



f3<p-f3- 1 (E ± -E)H v >J ± ±^ 



This implies 

HI (H - £±)<l>±||| 2 < e 2 + f3~ 2 (E± - Efe 2 



(H-E)ip 
f3~ 1 (E±—E)(H—E)ip 



C 2 e 2 . 



(3.3) 

(3.4) 
□ 



Remark 3.2. We know already from Proposition 2.10 (iii) that the spec- 



trum of H as given by (3.1 ) has always a gap between — \/3\ and \/3\ no 



matter how the spectrum of H looks like. For an interpretation of this 



in the context of Anderson's theorem, see Remark |2.11 (iii) 



We can also compute the density of states of H = (£ j 3 ^) , given the 
operator H possesses a density of states. The corresponding result is 
stated and proven here, even though we postpone the formal definition 
of the integrated density of states and of its Lebesgue derivative, the 
density of states, to the next section. 

Proposition 3.3. If H possesses an absolutely continuous integrated 
density of states with Lebesgue derivative D, then the integrated density 
of states of H is also absolutely continuous with Lebesgue derivative 
given by 

\E\ 



D(E) 



D{^W^W)+D{- JW^W) (3.5) 



VE 2 - /3 2 

for a.e. E e spec(H) C R\] - \(3\, andB(E) = for E £ spec(H). 
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Remark 3.4. Observe that in the case considered here the density of 
states D of H has a square-root singularity at the inner band edges ±/3 
whenever D(0+) 7^ 0, see the right panel of Figure [2] This is typically 
the case, if lies in the interior of the spectrum of a random Schrodinger 
operator H |Wl |Jj [HMJ . The square-root singularity may remind the 
reader of the van Hove singularity in dimension d — 1. In the present 
case the nature of the singularity is independent of the dimension d, 
however. 

Proof. By assumption we have for every interval An C R 

# (eigenvalues of lying in An} 
dE D(E ) = lim ^-^ — (3.6) 

A L ^°° |Al| 

where is a self-adjoint finite-volume restriction of H to the cube 
centred about the origin and containing |A^| = L d points for L G N 
odd. In particular, the spectrum of is discrete and converge 
to H in the limit L — > 00. Now, for a given interval A C [0, 00 [ we 
define 

A := {E eR: v^o + ^ e ^} (3-7) 

and 

H^:=( H(L) * V (3.8) 
V pi -hW J v 1 



Proposition |3.1| implies 

#{ eigenvalues of lying in An} 

= ^{eigenvalues of H < - L - ) lying in A}. (3.9) 
On the other hand, we have the identity 

/ dE D(E ) = fdE B(E), (3.10) 

which results from the change-of- variables E = \/ E 2 — f3 2 , E > 0, on 
A n [0, 00 [ and E = -^/E 2 - (3 2 , E > 0, on A n] - 00, 0]. Altogether 
we have 

, s #| eigenvalues of H^) lying in A} 

dE D(E) = lim — — f -. (3.11) 

a L ^°° |ArJ 



The same argument applies if A C] — 00, 0[. Thus, (3.11) generalises 



to arbitrary Borel sets AC1, proving the assertion. □ 

4. ERGODIC PROPERTIES OF RANDOM BLOCK OPERATORS 

^From now on we are concerned with random block operators. We 
consider a probability space (Q, J 7 , P) together with a block-operator- 
valued random variable 
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where FL is densely defined on H 2 = £ 2 {Z d ) © £ 2 {Z d ) for P-a.e. u e Q. 
Throughout we will assume that (at least) one of the situations of 



Proposition |2.1| applies to for P-a.e. w € 11 so that self-adjointness 
is ensured. 

We say that the random block operator H is ergodic (w.r.t. pL d - 
translations) , if there exists a period p G N d and an ergodic group of 
measure-preserving transformations {Tj}j £p z<i on Q such that H fulfils 
the covariance relation UjH^U* = H Tj ( w ) for every u G VL and for 

d / U \ 

every j G p1 d = X {Pk^), where Vj '■= [ n rr ) an d Uj is the unitary 
k=i \ j / 

translation operator on £ 2 (7* d ), that is Ujtp := <p(- — j) for every if G 

£ 2 {Z d ). 

Consequently, standard results [Kill ICLl \P¥\ IKi2j imply the exis- 
tence of a non-random closed set SCR, the a.s. spectrum of H, such 
that spec(H w ) = £ for P-a.e. u G Q. Analogous statements hold for 
the components in the Lebesgue decomposition of the spectrum. 

Next we introduce the central quantity of this paper which measures 
the density of spectral values of pZ d -ergodic random block operators 
H. 

Definition 4.1. The (non-random) right-continuous, non-decreasing 
function IN : R — >■ [0,1], defined by 

®(E) ■= ^7 E[tr^[XAoX]-oo^|(H)]] (4-2) 

for all E G R, is called integrated density of states of H. Here we in- 
troduced the elementary cell A := {j G Z d : ^ jk < Pk for all k = 
l,...,d}, E stands for the probabilistic expectation on Q, the trace 
extends over % 2 and the notation Xa refers to the multiplication op- 
erator ( ^ ^ ) where M is multiplication by the indicator function of 
A on £ 2 (Z d ). 

As we shall see from the next lemma, the definition of IN is justified 
by Birkhoff's ergodic theorem. For L G N odd, we denote by Al(j') the 
cube centred about j G 7L d and containing |Ai(j)| = L d many points 
of Z d . We also write A L := A L (0). 

Lemma 4.2. Let H be the random, F-a.s. self-adjoint block matrix 



operator (4.1), which is ergodic w.r.t. pZ d -translations . Then there 



exists a set VLq C Vt of full probability, P(f2 ) = 1, such that 

lim t 

L->oc 2 L d 



^( E ) = } [m ^fl tr « 2 [XArtM-aowOB*,)] (4.3) 



for every Bel, every u G an d every j 6 Z d . Moreover, the set of 
growth points offi is given by the a.s. spectrum ofM. 

Proof. This is fully analogous to |KM] or Sect. 5.1 in |Ki2j . □ 
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Mostly we will be interested in more specific random block operators 
whose diagonal blocks are given by the discrete random Schrodinger 
operator of the Anderson model in d dimensions and whose off-diagonal 
blocks are i.i.d. random multiplication operators. For simplicity and 
ease of presentation we dispense with the presence of magnetic fields. 
More precisely, let U : 7L d — > R be a periodic potential with period 
p EN d and define the bounded periodic background Hamiltonian Hq : = 
— A + Uo on H = £ 2 (Z d ). Consider two independent sequences of 
independent real-valued random variables (V(j)) . gZ<i and (b(j))j eZd 
on fl such that all the V(j) : oj !->■ V u (j) 's are identically distributed 
with law \iy and all the b(j) : oj b u (j) 's are identically distributed 
with law fib, and such that the V's are independent of the 6's. For 
the sake of technical simplicity we assume that the supports of the 
measures \iy and /i& are bounded subsets of E. For each elementary 
event oj G VL the functions V u : j i-> and : j i— >■ give rise 

to corresponding multiplication operators. Then the random operators 

H:uH>H u :=H + V u (4.4) 

and 

b : oj !->■ 6 W (4.5) 

are P-a.s. bounded and self-adjoint on £ 2 (Z d ). They are both ergodic 
w.r.t. pZ d -translations. From standard results, see e.g. |Kilt ICLt [PFj 
IKi2] . we know that their spectra are P-almost surely given by 

spec(H) = spec(H ) + supp(/iy) (4.6) 

and 

spec(6) = supp(/ife). (4.7) 
The corresponding random block operator 

H:oj^U w := ff _ b - ), (4.8) 



is also ergodic and P-a.s. bounded and self-adjoint on 7-L 2 . The next 
lemma roughly locates its a.s. spectrum. 



Lemma 4.3. Let H be the random block operator (4.8). Then we have 
P-a.s. the inclusions 

{ ± ^E 2 + (3 2 : E e spec(F),/3 G supp(// 6 )} C spec(H) C [-r,r] 

(4.9) 

where r := sup E&pec(H) \E\ + Bup^ 6Bupp(w) 

Remark 4.4. Since ^/ E 2 + (3 2 < E + (3, if both £, /3 > 0, we sus- 



pect from recalling Remark 2.11 (i) that even the left inclusion in (4.9) 



might be a strict one in certain cases. This is indeed true, as can be 
seen by choosing H = —A to be deterministic and //& a symmetrical 
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distribution w.r.t. reflection at zero. In this situation we may apply 

-a.s. 



Corollary |2.8| and get 

spec(H) = spec(H K ) = spec(A) +supp(/i5). 



(4.10) 



Combining Lemma 4.3 with Lemma 2.10 (iv) we infer 
Corollary 4.5. Assume that inf spec(H) ^ ; inf supp(/ib) ^ and 

A± := ±y/[in£ spec(H)} 2 + [inf supp(/i 6 )] 2 > 0. (4.11) 

Then A_ and X + are the endpoints of the open gap interval which sep- 
arates the positive and negative parts of the a.s. spectrum o/H. 



Proof of Lem ma\4.S\ While the right inclusion follows immediately 
from Lemma 2.10 (i) and (4.6), the left inclusion is based on a typi- 



cal Weyl-sequence argument. Let (y? n )neN be a Weyl sequence of nor- 
malised vectors for H and some fixed energy E G spec(ifo)- Since 
H is a bounded operator we may assume without loss of generality 
that each (p n ^ has compact support. The Borel-Cantelli lemma im- 
plies the existence of a set Qq of full probability, P(Oo) = 1> such that 
for every u G f2 , for every length L > 0, for every e > 0, for every 
v G supp Hv and every (3 G supp ^ there exist k G 7L d with 

\V u (j)-v\<e and \b u (j) - /3\ < e for all j G A^(fc). 

Now fix also v G supp /iy and (3 G supp We define 

E± := ±a/ (E + 

and, if /3 ^ 0, the vector 



+ /3 2 , 



(4.12) 



(4.13) 



± 



a/;. 



(4.14) 



where the constant M n ensures proper normalisation |||$^||| = 1. If 
(5 = 0, we set := ( v Q n ) and $~ := ( ° ). For every n G N we have 
the estimate 



|||(IL-£ ± )$ 



-{Vuj-v) 



Hq+v 

P 



P 

-(H +v) 



E. 



± 



■ n ' n 

and the identity 



(4.15) 



(4.16) 



(V u -v) 2 +(b u -/5 

For every n G N and every e > we exploit the freedom to shift the 
support of ip n such that it lies inside some large enough box A L (k) for 



which (4.12) holds. This implies 

, 2 



(Ti 1 ))^2, 2 ((^,$= 



2e 2 . 



(4.17) 
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On the other hand, given any e > 0, we infer from the Weyl-sequence 



property of (y? n )neN for H and Proposition 3.1 that there exists 



(2) 

n £ G N such that for all n ^ n £ the estimate T„ ^ Ce holds 
with some constant C > depending on E ,v and /3. This proves 
|||(H W — £? ± )$^||| < (C + a/2)£, and hence ($^) neN is a Weyl sequence 
for H and E±. □ 

In order to make manifest the interpretation of the integrated density 
of states IN as an eigenvalue counting function we have to introduce 
appropriate finite-volume restrictions of the block operator H. 



Definition 4.6. For H as in (4.8) and a finite cube Al C Z we 



introduce the Dirichlet and Neumann restrictions 

HSMf - ff V>) and HV:=C?K») (4.18) 



of H to the 2 1 A L | -dimensional Hilbert space H\ := £ 2 (A L ) © i 2 {A L ). 
We also introduce the Dirichlet-bracketing and Neumann-bracketing 
restrictions 



Hr-r"-^J and H - : =rr-<o- ( 4 - 19 ) 

In the above we have used the Dirichlet, resp. Neumann, restriction 

H%) N := -A^ + U + V (4.20) 

of H to £ 2 (Al), and all multiplication operators are to be interpreted 
as canonical restrictions. 

Remarks 4.7. (i) The Neumann Laplacian on £ 2 (Al) is also called 
graph Laplacian or combinatorial Laplacian. It is defined by 

(-A£ty)(j):= ^ ( 4 - 21 ) 

fceAz,:|j-fc|=l 

for every -0 G £ 2 (A^) and every j G A^. The Dirichlet Laplacian can 
be represented as the perturbation 

- := -A^ } + 2r (L) , (4.22) 

where r^ L ^ is the multiplication operator given by 

(I^VXj) := VC?) G Z d \ A L : \k - j\ = 1} (4.23) 

for every j E A L . In other words, lives only on the outermost layer 
of Al and multiplies by the number of missing neighbours at every 
point. 

(ii) The finite-volume restrictions — A^ and — A^ are defined 
such as to obey Dirichlet-Neumann bracketing. Consequently we obtain 
the chains of inequalities 

hL L) ^ and hL L) < ^ (4.24) 
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for every Al C Z d and the Dirichlet-Neumann bracketing properties 

H (A) e H (A') ^ H (AUA') H (A) e jjCA' ) ^ jjCAuA') ^ 

for all disjoint A, A' C Z d (in obvious abuse of our notation). 
The desired interpretation of IN follows from 



Lemma 4.8. Given the random block operator H from (4.8) and any 

E G R ; we define the random, finite-volume eigenvalue counting func- 
tion 

W ( i\E) := tr«. [XHoo^CH^)], (4.26) 

where X symbolises any self-adjoint restriction such that 

M {L) <: < H^ L) . (4.27) 

T/ien 

(i) there is a set Q C Q of full probability, P(fio) = 1; suc/i 

1N(£) = lim N<g,(£) (4.28) 
/or even/ cj e f2 0; even/ continuity point E e 1R o/lN and ever?/ bound- 



ary condition X satisfying (4.27). 
(ii) In addition we have the bounds 

E[]N^ L) (£)] < M(£) < E[M^ (E)] (4.29) 
/or every finite cube A^ and every 

Proof. The arguments proceed as in the standard random Schrodinger 
case, see e.g. [HH ED EH EMI E2] • □ 

5. Random diagonal B: boundedness of the density of 

STATES 

We have seen in Section [3] that the density of states of H with off- 
diagonal operators B that are a constant multiple of the identity can 
have a singularity, even if the density of states of H is bounded. In this 



section we consider random operators H of the form (4.8) with diag- 
onal disorder in all blocks. We will prove a Wegner estimate for such 
models. It implies that the density of states of H exists and is bounded 
provided H or b is bounded away from zero and the distribution of 
the random variables is absolutely continuous with a suitably regular 
Lebesgue density. Technically, one has to cope with a non-monotone 
dependence on the random potential: the single-site potential that en- 
ters H is sign-indefinite and has mean zero. In contrast, for ordinary 
(i.e. non-block) random Schrodinger operators with such single-site po- 
tentials it is not known how to prove a comparable Wegner estimate. 
Existing Wegner estimates involve either higher powers of the volume 
|Ve] or yield only log-Holder continuity of the integrated density of 
states [El!. 
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Figure 3. Assuming b u (j) = l+wb u (j) with b u (j) dis- 
tributed uniformly in the interval [—0.5, 0.5] we display 
the dependence of the density of states on the disorder 
parameter w. (Assumption (B) of Theorem 5.1 is satis- 
fied for w = 0.1,0.5,0.9). 



Theorem 5.1. Consider the random block operator 

H : u ^ R, = 



(5.1) 



where H and b are given as in (4.4) and (4.5). Assume that at least 



one of the following two conditions is satisfied: 

(H) there exists A > such that H ^ XI holds F-a.s. and fly is ab- 
solutely continuous with a piecewise continuous Lebesgue density 
cfiv of bounded variation and compact support, 

(B) there exists (3 > such that b ^ (31 holds F-a.s. and fib is ab- 
solutely continuous with a piecewise continuous Lebesgue density 
<fib of bounded variation and compact support. 

Then the integrated density of states IN of H is Lipschitz continuous 
and has a bounded density D := dM/dE. In particular, we have for 
Lebesgue- almost all E e IR 



D(E) *C 2 



\E\ + 1 
A 



V||bv> 



(5.2) 
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if Condition (H) applies, and 



D (£) <2 ^L^||0 6 || BV) (5.3) 



if Condition (B) applies. Here, ||/|| B v denotes the total variation of 
f : K -»■ C. 

Remarks 5.2. (i) The density of states D is an even function on 
account of Lemma I 



(ii) We recall that, typically, D exhibits a singularity at the inner 
band edges for constant off-diagonal blocks b — /31 > 0, see the left 
panel of Figure [I] Case (B) of Theorem 5.1 implies that, regardless 



of iJ, this singularity is smeared out by the disorder. In particular, D 
remains bounded. 

(iii) Figure [3j displays a typical density of states for random block 



operators in the case (H) of Theorem 5.1 Here, D is not only bounded 
but even vanishes at the inner band edges, because it exhibits a Lifshits 
tail, see Section |6l 



Proof of Theorem 5.1, As compared to the proof of the Wegner esti- 
mate for the standard Anderson model, there are essentially two mod- 
ifications necessary here. They have been isolated in the subsequent 
Lemmas 15.31 and 15.41 



Lemma 2.2 tells us that the roles of the diagonal operator H and 



of the off-diagonal operator b are interchangeable. Thus it suffices 



to prove only the case (H). Moreover, because of Lemma 2.3 we will 



restrict ourselves to energies E ^ without loss of generality. In fact, 



using Condition (H) and Proposition 2.10 (ii) it then suffices to consider 
E ^ A, which we will do from now on. We denote by H( L )(V) : = 
the Dirichlet restriction ( |4.18[ ) of the operator H to the cube A^. 
In this notation we make explicit the dependence of the operator on 
the random variables V := {Vj}j(z\ L (we prefer to write Vj instead 
of V(j) in this proof and the subsequent lemmas) and suppress the 
dependence on the entries of the off-diagonal blocks. Furthermore, we 
write E n (V) := E n (M^(y)) for the n th eigenvalue of H( L )(V), where 
the eigenvalues are ordered by magnitude and repeated according to 
multiplicity. Finally, we fix e G]0, min{A, l}/3[ and consider a switch 
function g, i.e. g G C 1 (M) is non-decreasing with ^ g ^ 1, g(rj) = 1 
for 1] > e and g(rj) = for rj < —e. Then 



< X]E-e,E+ £ ](v) < q(v ~ E + 2e) - q{t] - E - 2s) 



(5.4) 
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for all 77 e R, whence 

W [x ]S - £lS+£ ](H( L )(V)) 



2|Ai| 

^ ^ [g(E n (V) - E + 2e) - f?(£ n (V) — E — 2e) 

n=l 

2|Ai| £ +2e 



E 



5 



2|A L | s+2£ 



e(s n (V) - 77) d?? 



E 



e'(S n (V)-77) dry. 



n=l JS " 2£ 

The chain rule tells us that 



(5.5) 



(5.6) 



jeA L J jeA L 
Since H^{y) ^ A > by Condition (H), we conclude from this iden- 



tity and Lemma 5.3 



A 

E + l 



9 



J6A L J 



(5.7) 



for all n G N and all 77 e [J5 - 2e, £ + 2e]. We note that the last 
inequality uses 3e < min{A, 1}. This guarantees that only those n 
with E n (V) e]0, E + 1[ contribute and that the j-sum is non-negative 



for these n by Lemma 5.3 



Since the random variables {Vj}j & \ L are independent and have the 
same individual distribution fiy, the expectation is just integration 
with respect to the product of these distributions. Thus (5.5) and 
(5.7) imply 

E{tr^2 [X]E-e,E+e](H 
E + l 



(L) 
D 



A 



E 



E+2s 
E-2e 



I r a 2|Ail 

x ( 11 d M^fc)W I 



keA L : k+j 



Since, in general, the function 

2|Ai| 

y. _> F{ y.) := g(E n (V) -rj)= t %| Uh^(V) - V ) (5.9) 



n=l 
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is non-monotone in its argument for given r] G R and Vk G R, 7^ j, we 
deviate from the standard reasoning at this point. Clearly F G C 1 (1R) 
by analytic perturbation theory. Moreover changing Vj amounts to a 
rank-2-perturbation of H^ L ^(V). Thus we have \F(v) — F(v')\ ^ 2 for 
all v,v' G R and we can apply Lemma 5.4 to (5.8). This gives 

E{tr wi [XiE-cE+e]^)]} < 8eL d \\ ( j )v \\ BV E ^. (5.10) 

Finally, the assertion follows from Lemma 4.28 and dominated conver- 
gence. □ 

We come to the main deterministic tool used in the proof of Theorem 



5.1 Given Assumption (H), it ensures that eigenvalues move around 



strongly enough when the random variables change their values. 
Lemma 5.3. Let E(V) be an eigenvalue of ^ b{l) H(L) ^j. Then 

E(V) > inf s P ec (^ (L) (V)). (5.11) 

Proof. Let (V) = = (£) G U\ be an eigenvector corresponding to 
the eigenvalue E(V) = E, normalised according to ((^, ^/)) = (01, ipx) + 
(^2,^2) = !■ Thus > writing ffW(V) = it satisfies 

# (L Vi + 5 (L V 2 = E^x, 

5 (L Vi - # (L V 2 = Eip 2 . 

By the Feynman-Hellmann theorem (see e.g. [RS]) we have 

dE 



(5.12) 



dV 3 



\MJ)\ 2 ~ \MJ)\ 2 (5-13) 



for every j G Al and consequently, inserting (5.12) we find 

= (Vi, # (L Vi + £ {L V 2 > - (^ (L Vi - # (i V 2 , 

^ inf spec (tf (L) ). (5.14) 

□ 

The next lemma deals with the problem of the non-monotonous de- 
pendence of the cumulative eigenvalue counting function on the random 
potential. It is here where we have to assume a suitable regularity of 
the distribution. 

Lemma 5.4. Let : R — > C be piecewise continuous, of bounded 
variation and have compact support. Let F G C 1 (IR) and assume the 
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existence of a constant a 
i,!/Gl. Then we have 



> such that \F(x) — F(y)\ ^ a for all 



F'(x) 4>(x) dx 



^ a 



(5.15) 



Proof. Step 1: We prove the claim for step functions 

N N 

0JV = ^2^uX] Xv -i^] = ~ ^~^W-i,x N ] (5-16) 

v=\ u=l 

where N G N, £ v G C for all v G {1, . . . , iV}, £ := and -oo < x < 
x\ < . . . < xn < oo. For such 0^ we get 

N 



F'(x) 



'x) dx 



^(e,-e,-i)[^iv)-F(x,- a )] 



^=1 

A? 



a \\<p N \ 



(5.17) 



i/=i 



Step 2: For general as in the lemma, the claim (5.15) follows from 
a uniform approximation of by step functions of the form 0jy := 
Siv=i 4 > { x v)'X.]xu-i,xu\- Indeed, for every given e > we can choose 
discretisation points x u , v = 1,...,N and iV G N such that ||0 — 
<Pn\\oo ^ e, because is piecewise uniformly continuous. Therefore 
L F'(x) <f>(x) dx = limAr.j.oo J R F'(x) 4>n(x) dx for a suitable sequence 
of step functions 0at, because F' is bounded on the support of 0. This, 
the bound from Step 1 and ||0jv||bv ^ II^IIbv finish the proof. □ 



6. Random diagonal B: Lifshits tails 

It is a striking fact that the integrated density of states of random 
Schrodinger operators grows only exponentially slowly in the vicinity 
of fluctuation band edges. This behaviour is called Lifshits tail [L]. 



In Theorem 6.1 we provide a result for random block operators H 
of the form (4.8) which limits the growth of the integrated density of 



states IN at energy A, provided H ^ Al > exhibits a Lifshits tail at 



energy A. We emphasise that Theorem 6.1 is only interesting in the case 
where ±A coincide with the endpoints of the spectral gap of H around 



zero. According to Corollary 4.5 this always happens if G supp (//&). 



Theorem 6.1. Consider the random block operator 

H : co ^ H w = 



b u -H u 



(6.1) 



where H and b are given as in (4.4) and (4.5) and let IN be its integrated 
density of states. Suppose in addition that A := inf spec(if) > and 
that there exists constants a, 7 > such that for all sufficiently small 
s>0 



P[inf spec(#<^ e) ) < A + e] < e" 



(6.2) 
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where L £ is a sequence of diverging lengths as e I with lim e ^ £ a ^ d L £ 
exists and lies in ]0, oo[. Then the estimate 



lim sup 

e4,o 



In | ln[F(A + e) - 1N(A)] | 



sC -a 



(6.3) 



holds. 



Remarks 6.2. 



Assumption (6.2) is the statement which is typ- 



ically proven for a random Schrddinger operators H of the form (4.4) 
when establishing the upper bound 



In I lnLY(A + e)] I 

lim sup — : ^ 

£ ;o In £ 



-a: 



(6.4) 



for a Lifshits tail with Lifshits exponent a = d/2 of the integrated 
density of states N at the lower spectral edge A. In other words, if 
G supp(/ib), then Theorem 6.1 says that the growth of IN near the 
lower edge of the positive a.s. spectrum of H is no faster than the 
growth of N near the bottom of the a.s. spectrum of H. Note that the 
bound on IN holds independently of b. 



(ii) An analogous statement to Theorem 6.1 holds at the upper 
edge of the negative spectrum of H. 



Proof of Theorem (LI . We claim that for all energies E ^ and all 
finite cubes A L the estimate 



N(E) - 1N(0) ^ E[1N W (£) - ]NL L, (0)] 



(6.5) 



holds. This follows from Dirichlet-Neumann bracketing (4.29) and 

1 

2 



1N(0) 



E[1N 



A' 



EpNl L) (0)]. 



(6.6) 



The first two equalities in (6.6) are based on Lemma 2.3 and that zero 
is not in the spectrum. In order to see the last equality in ( 6.6[ ) we view 



H 



H( L )(1) as an analytic perturbation of = H^^O), where 



HW(a) :=U { N L) -2a(° o 



o o 

r(i) 



(6.7) 



for a e R and was introduced in (14.221). Analytic perturbation 
theory tells us that the eigenvalues of H^(a) depend continuously 
on the parameter a. On the other hand, we infer from Lemma |2.12| 
that zero lies in an open spectral gap of TBS L '(a) of size at least 2A for 
every a E [0, 1]. Therefore H^(0) must have exactly as many positive 
(negative) eigenvalues as H^ L ^(1), and the third equality in (6.6) holds. 
Now let e > and observe 



E[lNL i) (A + £) -^(0)] < P[spec(H (Lj )n]0,A + e] ^ 0]. (6.8 
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By contradiction we conclude from Lemma 2.10 (ii)| that the event in the 
probability on the r.h.s. of (6.8) implies the event inf spec(if^) ^ X+e. 



Hence, we get 

1N(A + e) - N(0) < P [ inf spec (H^) ^ X + e], (6.9) 



and the claim follows from (6.2). □ 
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